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Abstract 



| Nonlinear perturbations of Friedmann-Lemaitre cosmologies with dust and 

• a cosmological constant A > have recently attracted considerable attention. 

. In this paper our first goal is to compare the evolution of the first and second 

' order perturbations by determining their asymptotic behaviour at late times 

in ever-expanding models. We show that in the presence of spatial curvature 
K or a cosmological constant, the density perturbation approaches a finite 
^ , limit both to first and second order, but the rate of approach depends on 

the model, being power law in the scale factor if A > but logarithmic 
if A = and K < 0. Scalar perturbations in general contain a growing 
and a decaying mode. We find, somewhat surprisingly, that if A > the 
decaying mode does not die away, i.e. it contributes on an equal footing as 
the growing mode to the asymptotic expression for the density perturbation. 
On the other hand, the future asymptotic regime of the Einstein-de Sitter 
universe {K = A = 0) is completely different, as exemplified by the density 
perturbation which diverges; moreover, the second order perturbation diverges 
faster than the first order perturbation, which suggests that the Einstein- 
de Sitter universe is unstable to perturbations, and that the perturbation 
series do not converge towards the future. We conclude that the presence of 
spatial curvature or a cosmological constant stabilizes the perturbations. Our 
second goal is to derive an explicit expression for the second order density 
perturbation that can be used to study the effects of including a cosmological 
constant and spatial curvature. 
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1 Introduction 

In a recent paper Uggla and Wainwright (2013) (referred to as UW3)0we gave a sim- 
plified structure for the system of equations that govern second order cosmological 
perturbations. Because of their generahty these equations provide a starting point 
for determining the behaviour of nonhnear perturbations of Friedmann-Lemaitre 
(FL) cosmologies with any given stress-energy content, using either the Poisson 
gauge or the uniform curvature gauge. In the present paper we use this system of 
equations, specialized to the Poisson gauge, to analyze the behaviour of linear and 
nonlinear perturbations of FL cosmologies with dust and possibly a cosmological 
constant as the matter-energy content. This problem has a lengthy history starting 
with Tomita (1967), and has received considerable attention over the past fifteen 
years0 

Our first goal is to compare the evolution in time of the nonlinear (second order) 
perturbations to the linear (first order) perturbations by determining the asymptotic 
behaviour of the perturbations at late times in ever-expanding models, which has 
received little attentionjf] Since we do not wish to make the customary assumption 
that the spatial geometry of the background is flat 3 there are three cases to consider: 

i) A > 0,K arbitrary, which we refer to as the de Sitter asymptotic regime, 

ii) A = 0, K < 0, which we refer to as the Milne asymptotic regime, 

iii) A = 0, K = 0, which we refer to as the Einstein-de Sitter asymptotic regime. 

We show that in the de Sitter and Milne asymptotic regimes, the second order 
perturbations have the same asymptotic time dependence as the first order pertur- 
bations. For example, in both cases the density perturbation "freezes in", i.e. it 
approaches a finite limit as x — t- oo, where x = a/oo is the dimensionless background 
scale factor: 

hm = (^•)(5oo, (1) 

a;— >oo 

^The present paper also relies heavily on Uggla and Wainwright (2011) and Uggla and Wain- 
wright (2012), which we shall refer to as UWl and UW2, respectively. 

^See for example, Matarrese et al (1998), Bartolo et al (2005), Tomita (2005), Boubekeur et al 
(2008), Bartolo et al (2010) and Hwang et al (2012). 

■^The only papers of which we are aware are Bruni et al (2002) and Mena et al (2002). They 
consider non-linear perturbations of a flat FL universe with dust and a positive cosmological 
constant. Their approach is heuristic in that they solve a truncated asymptotic version of the 
evolution equation for the perturbations in the synchronous gauge, whereas we solve the general 
evolution equation in the Poisson gauge. In addition we focus on the matter density perturbation 
at second order whereas they deal with the metric perturbations. 

''Current observations restrict Vtk to be close to zero. See, for example, Okouma et al (2013). 
However, since a flat FL model is unstable to perturbations of VL^ away from zero it is unlikely 
that Vtk will remain close to zero. Thus we are interested in how the presence of spatial curvature 
affects the future asymptotic behaviour of perturbations. 
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where r = 1,2 labels the first order and second order perturbations, respectively. 
There are two significant differences, however. First, in the de Sitter regime, the 
rate of approach to the asymptotic state is power law in x, while in the Milne regime 
it is logarithmic: 



de Sitter regime it depends on both the growing and decaying modes of the linear 
perturbation, while in the Milne regime it depends only on the growing mode at 
first order. In other words, in the de Sitter asymptotic regime the decaying mode 
does not die away. 

Finally, we show that the Einstein-de Sitter asymptotic regime is totally different 
from the other two regimes. This is typified by the density perturbation, which 
diverges: 



In particular, the unbounded growth of these physical perturbations in the Einstein- 
de Sitter asymptotic regime suggests that the Einstein-de Sitter cosmology is un- 
stable to perturbations, and that the perturbation series will not converge for x 
sufficiently large. On the other hand our analysis of the de Sitter and Milne asymp- 
totic regimes as described above, shows that the presence of spatial curvature or of 
a cosmological constant stabilizes the perturbations. 

Our second goal is to derive an explicit expression for the second order density 
perturbation that can be used to study the effects of including a cosmological con- 
stant and spatial curvature. In order to obtain a relatively simple expression which 
nevertheless illustrates various interesting phenomena, we consider the special case 
in which the decreasing mode of the linear perturbation is assumed to be zero, as 
is usually done in cosmological perturbation theory. When the spatial curvature is 
set to zero our expression is related to that given recently by Bartolo et al (2010). 

The outline of the paper is as follows. In section 2 we give a new and concise 
derivation of the solutions of the equations that govern first order and second order 
perturbations, expressing them in integral form. Our derivation is made simple by 
our use of the factorization property of the linear differential operator that governs 
the evolutionjf] We show that if the decaying mode of the scalar perturbation is 
set to zero the general solution can be given in explicit form, which leads to the 
expression for the density perturbation referred to above. In section 3 we determine 
the asymptotic behaviour of the perturbations in the three asymptotic regimes, and 
draw the conclusions described above. Section 4 contains a brief summary and 
discussion. The details of the analysis of the source terms and of the asymptotic 
behaviour are given in the appendix. 

2 Perturbed dust cosmologies to second order 

In this section we solve the governing equations for linear and second order per- 
turbations of dust cosmologies, using the Poisson gauge. We make the simplifying 





as X — oo. 



(3) 



^See UWl equation (55) and UW2 equation (41). 
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assumption that the perturbations at the linear order are purely scalar. This imphes 
that the matter has zero vorticityj^ 

2.1 Background 

We consider perturbations of a FL cosmology containing pressure-free matter (dust) 
and a cosmological constant A > 0. The background Robertson- Walker metric is 
given by 

ds"^ = a^{-drj^ + jijdx'dx^), (4) 

where t] is conformal time, a is the scale factor and 'jij a 3-metric of constant 
curvature with curvature index K. The scale factor a determines the dimensionless 
background Hubble scalar "H according to 

U = - = aH, (5) 

a 

where H is the true background Hubble scalar and a prime denotes a derivative with 
respect to rj. As in UW3 we introduce the scalars 

AG = 2{n'-H^ + K), ^^ = a2((o)p + (o)j9), (6) 

where 

(°)p = (V + A, (°V = ^°W-A, (7) 
where the subscript m refers to the matter component. We assume that the Einstein 
equations in the background are satisfied, which imply that Ac = At- We thus drop 
the subscripts G and T, and for pressure-free matter we have 

A = 2{H' -'H^ + K) = a^ (V- (8) 

We introduce a reference time 7]q and a normalized scale factor x defined by 

X := a/ao, where = a(?7o)- (9) 

For a dust source the matter conservation equation in the background yields (x^^^^p)' = 
0, which implies that (Ax)' = 0. We define a constant m by 

\Ax. (10) 



,2 



3 



We also introduce the usual density parameters: 



(°)pm K A 

:= '^jp^ '■= •= '^Jp' ^^^^ 

and note that for dust 

A = 3'H''nm, m^ = 'Hlnm,o. (12) 

The Friedmann equation in the background, Vtm + fifc + f^A = 1, can now be used 
to express "H explicitly function of x: 

= nl (fiA.ox^ + ^kfl + n„,,ox-^) , (13) 

with the subscript indicating evaluation at 7]q. 



'See equation (B.41d) in UWl. 
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2.2 The governing equations 

In the Poisson gauge the gauge invariants that describe the metric perturbation are 
the Bardeen curvature and the Bardeen potential and those that describe 
the matter perturbation in the case of pressure- free matter (dust) and a cosmological 
constant are ^^^S, the density perturbation, and '■''V, the velocity perturbation, where 
r = 1 at first order and r = 2 at second order. We refer to UW3 for the definition of 
and (see equation (25)), and to equations fll02p and (11041) in the appendix of 
this paper for the definitions of ^"^^d and *^^V. In UW3, in order to write the governing 
equations for first and second order perturbations in a so-called minimal form we 
also defined matter gauge invariants in terms of the components of the stress-energy 
tensor, denoted by '■^-ID) and '•''V, which we give in equations (I106P in the appendix. 

We now specialize the governing equations derived in UW3l3 to the case of dust. 
The Bardeen curvature plays a central role in determining the dynamics of the 
perturbations. At each order it satisfies a second order linear differential equation, 
which we write in operator notation in the following form: 

= 0, = S, (14a) 

where the differential operator C is given by 

C = dl + 2,nd^ + {2n' + V? -K), (14b) 

and the source term S depends quadratically on and its derivatives. We thus 
have to solve these equations successively, with the solution of the first equation 
determining the source term on the right side of the second equation. The remaining 
perturbation variables ^^to, "^^'V and with r = 1, 2, are determined by and 
^^KE' either algebraically or by differentiation, as follows: 

(% = 2^-i(D2 + 3ir)(%, (% = 2^-i(D2 + 3ir)(% + 5o, (15a) 

H(% = -2^-^H(5^ + H)(%, H^^^ = -2A-^H{d^ + Hi^^ + Sy, (15b) 
= (^)$ = (2)^ + 5$, (15c) 

where = 7*''DjDj and Dj is the spatial covariant derivative associated with 
7jj. The complete expressions for the source terms iS,iSid),iSv and 5$, which depend 
quadratically on and its derivatives, are given in equations fl70p in the appendix. 
Once (''to and (''V have been calculated one can obtain '^'^^5 and '■^V by using the 
following relations, that are a special case of equations f lllOp in the appendix: 

(1)5 = (% + 3?/(%, ^% = + m^^^ - 2(D(%)2, (16a) 

= (2)^ = - 25* [((i)(5 - (%)D,((i)v)] . (16b) 

For our purposes it is important that the operator C in equations ffT^ can be 
written as the product of two first order differential operators: 

£(.) = HCaCb (^) , (17a) 



^See equations (61). 
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where 

-HCa := ndr, + 2^.' + CB:=d^ + 2'H, (17b) 

(see UWl, equation (55) and UW2, equation (39)). Replacing conformal time r] by 
the normahzed scale factor x, noting that d„ = Tixdr, we can now write C in the 
torno 



2 

X 



2.3 Linear perturbations 

The first order Bardeen curvature ^^^"^ is obtained by solving the first of equa- 
tions fll4ap . = 0. We can use ffTS]) to integrate this equation twice, obtaining 
the general solution in the formj^ 

^ix,x') = ^ {C+{x')Iix)+C.{x')) , (19a) 



x^ 



where 



■= I (19b) 



'H(x 



The spatial functions C+(x*) and C_(x*) describe the growing and decaying 
modes of the perturbation, respectively. The function C+(x*) is of particular sig- 
nificance, for the following reason. Two "conserved quantities" that are associated 
with scalar perturbations of FL have been defined in the literature. These quanti- 
ties, often denoted by the symbol ( with a subscript, satisfy an evolution equation 
of the form 

a^C. = D2C., (20) 

where C, is an expression involving the primary gauge invariants such as \I' or 
V and the background variables. This equation is referred to as a "conservation 
law", since if spatial derivatives are negligible ("perturbations outside the horizon") 
in some epoch, then fl20l) is approximated by dri(, = 0, i.e. (, is approximately 
constant in time during that epoch. We refer to UW2, section 4, for a unified 
discussion of these conserved quantities. In particular the conserved quantity that 
we denote by Cv can be expressed in terms of ^ and V according tc0 

Cv = f 1 - ^1 ^ - nY. (21) 



A J 

On substituting for V from the first of equations fll5bp we can rearrange fl2T]) to 
obtain 

Cv = ( ■ (22) 



Ax Kn 



^Write the operators Ca^Cb in the form HLAi*) = dxCH^x*) and >Cs(») = ^dx{x^ •)■ 
^We will henceforth drop the superscript on first order quantities whenever there is no danger 
of confusion. 

^"^See UW2, equation (72). This quantity is sometimes referred to as "the curvature perturbation 
in the comoving gauge". See Malik and Wands (2009), equation (7.46), who use the symbol TZ for 

Cv. 
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Substituting for \1/ from fll9ap and using (fTU]) leads to 

Cv = lm-'C+. (23) 

This result shows firstly that for any perturbation with dust, the conserved quantity 
Cv is exactly constant in time, and secondly it relates the arbitrary spatial function 
C+(x*) to Cv We can thus write the general solution fll9al) in the form 

^{x,x') = fm^^ (/(x)C(xO + C_(xO) , (24) 

on rescaling C_(x*). Here and in the rest of the paper, for notational convenience 
we drop the subscript v on We note in passing that in other discussions of 
scalar perturbations of dust, another conserved quantity, the so-called curvature 
perturbation on constant density hypersurfacesl"1 plays an important role (see, for 
example, Bartolo et al (2006), equation (2.13)). 

The other gauge invariants can be expressed in terms of First, the matter 
gauge invariant D is given by the first of equations ( 115ap . Second, the gauge invariant 
V is given by the first of equations ( llSbj) . or can be calculated using ( 12T1) since C, is 
constant in time. Finally we have $ = \& as follows from the first of equations fllScI) . 



2.4 Second order perturbations 

The second order Bardeen curvature is obtained by solving the second of equa- 
tions fll4ap . = iS(\l/). The solution of this equation that satisfies the initial 
conditions 

(%(x,„,i,xO = 0, 9,(%(x,„,i,xO = 0, (25) 

where Xinu > 0, can be obtained by using f|T8|) and integrating twice. This leads to 
the following formula: 

(%(x,x^) = ^ r ^^dx, (26a) 

where we have defined ^ 

S(x,x^) := / S{x,x')dx. (26b) 

Here 

S{x,x') = S{^{x,x')), (26c) 
with \I^(a;,x*) given by (!24l) . The general solution can be written in the form 

(%5en = + ^%hon^o„ (27) 

where ^'^^ homog is the general solution of the homogeneous equation = and 

hence is of the form fll9ap . with C+(a;*) and C_(x*) being determined by the values 
of ^"^^gen and drj^'^^gen at the initial time Xmu- 

^^This quantity is denoted by ( in Malik and Wands (2009) (see equation (7.61)), and by Cp in 
UW2 (see equations (65)-(67)). In the long wavelength regime Cp and Cv are approximately equal, 
but unlike ^v, Cp is not exactly constant at first order. 
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The other perturbed quantities at second order, '^^-'D, and are determined 
exphcitly in terms of ^"^^ and by the constraint equations (fT5|) . We now give 
an alternative expression for ^^V. By differentiating fl26ap with respect to r], using 
drj = Hxdx, one can express (9^ + "H)*^^^ in terms of E>{x,x^) and This leads 
to the following expression 

n = -|m-2 §(a;, x') + (^1 - (% + 5v, (28) 

which is useful for analyzing the asymptotic behaviour of ^^V. 

Finally we note that the solution f l26al) can be written in an alternate form by 
changing the order of integration, yielding 



(%(x,x*) = ^ / [/(x) - I{x)]Six,x')dx. (29) 



2.5 Zero decaying mode at linear order 

The solution given by fl26al) or (!29|) . depends on the two arbitrary spatial 
functions ( and C_ through the source term S which depends on We can obtain 
detailed information about the time and spatial dependence of '•^Hl/ by considering 
the special case of a linear perturbation with zero decaying mode^ i.e. C_(x*) = 0. 

When the decaying mode is zero, equation IHM and the first of equations (115bll 
reduce to 

vl/(x, x') = g{x)C{x^), = -|fi„V^7C, (30) 



C, (31) 



wheri 



/(a:):=l + ^, g{x):=lm'%. (32) 
g x^ 

We use (1291) to obtain a particular solution for In this case one can if desired 
choose Xinit = in fl29|) . since S = 0(1) as a; — )■ oo, as follows from fl7^ and fl75l) . 
We now substitute the expression (1741) for the source term 5 into fl29|) . which leads 
to 

(%(a;, = ^ (5i(x)C' + B2{x)V{0 + m-^[B;{x){BCy + B4x)D^V{C)]) , 

(33a) 



^^This simplifying assumption is often made, usually because it is claimed that the decaying 
mode will become negligible compared to the growing mode in the future. However, as we show 
in Section rS.l.ll this is not the case if A > 0: the decaying mode contributes to the perturbations 
at both linear and second order. If one assumes that the decaying mode is zero one is essen- 
tially considering perturbations in a universe with an isotropic singularity (Goode and Wainwright 
(1985)). 

^■^We note in passing that the function := xg{x), defined up to a constant factor, is sometimes 
referred to as the growth suppression factor. See for example Bartolo et al (2006), in the text 
following equation (2.3). In our analysis the related function I{x) plays a central role. 
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where 



Ba(x) 



n 

X 



[I{x)-I{x)]TA{x)dx, A = l,...,4, 

-- init 



(33b) 



and the functions Ta{x) are given by (175|) . The notation (D/)^ and the operator V 
are defined in (1721) . 

The density perturbation ^'^^S is given by fll6ap . We calculate the expression 
+ SH^"^^ that is required by substituting as given by fl33|) . into the second 
of equations ( I15ap and (I15bp . After using f lTOj) . f|T2|) and 9^ = "Hx^^,. we obtain 



+ 3-H(% ^ _2^];;^l4(x(%) + |m-2D2(a;(%) + 5d + 35v. 
The source terms So and iSy are given by (1701) . The final result is 



where 





-2n 


A2 = 


-2VL 


A3 = 


-2VL 


A,= 


-2VL 



4.5i-/((l-/)2-4fifc) 

4■52-V(l + |^]-V')" , 



-1 



with 



A., = \{B^ + Axg^ 



(34) 

(35) 

(36a) 
(36b) 
(36c) 
(36d) 

(a + ^^fcX-^)5A. (36e) 

Here the time-dependent functions Ai and A2 identify the Newtonian terms that 
dominate at late times, ^3 — A^ identify the post-Newtonian terms, while S3 and 
B4 identify the super-horizon terms that describe the perturbations on the largest 
scales. These time-dependent functions depend on A and K through the functions 
T-Li^x) and I{x), and when K = = A and we choose Xmu = they reduce to powers 
of X, as given in section 12.5.11 Observe that the spatial dependence is described by 
expressions such as (DCY that are quadratic in ( and its spatial derivatives. In the 
super-horizon, post-Newtonian and Newtonian terms these expressions are of degree 
zero, two and four in the spatial derivative operator D, respectively. 

Equation f p3|) gives a particular solution for (2)v]> that satisfies lim^-^o ^"^^ = 0, 
and fl35p gives the corresponding expression for the density perturbation ^'^^d. The 
general solution for ^^K]/ subject to the condition that the decaying mode at linear 
order is zero is given b\l^ 



gen 



(% + C(xO^(x) 



(37) 



where ^^Kl' is given by ( l33l) . The corresponding density perturbation is given by 



(2)^ 



gen 



+ Is 



X ■ 



C-31]~i(/ + fi,)C 



(3^ 



*The term (2)^ 



honiog 



in (HZl) is of the form (ITOal) with C_ = 
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where ^"^^S is given by (135|) . 

The expression for ^'^'^dgen given by (l35l) . (136|) and (l38l) is new and represents 
one of the main results of this paper. We have shown that the expression for ^'^^d 
given by Bartolo et al (2010) when K = and A > can be written in our form0 
provided that the spatial function C{x^) is chosen suitably. Specifically, in deriving 
their result they use the level of primordial non-Gaussianity at the end of inflation 
{ibid, section 3) to determine the function C(a;*) in f l37|) . which has the form 

C{x') = lgr42V{0 + (1 - fanOCl, (39) 

where a^i is a constant that parametrizes the primordial non-Gaussianity and 
is the value of g at some initial time. Finally we note that Tomita (2005) has also 
given an expression for ^"^^S when K = and A > 0, in which the time dependence 
is expressed in a completely different way@ 

2.5.1 A special case: Einstein-de Sitter 



When = and A = 0, it follows from ([TS]) and (fTObj) that "H and / {x) are given 
by 

= m^x-\ mH{x) = lx'''\ (40) 
and hence that the time dependence functions / and g are 

g{x) = I, fix) = 1. (41) 

One can now use ( 1751) to calculate the functions Ta{x), which when substituted 
in ( I33b|) with Xmit = leads to 



fii = = B2, B; = ^^x\ 54 = 2053. (42) 
Equations (!33|) . (!35l) and (136!) now lead to 

(% = ^(C)x, with ^(C):= if5m-2[(DC)2 + 20D2p(C)], (43) 

i% = Mp(^) + _|_^-2 (137(DC)' - 80D2i?(C) + 840^(2) x+lm-^T)^F{C)x^ . (44) 

The corresponding general solution is given by ( 1371) and ( l38l) . The second order 
density perturbation in the Poisson gauge for the Einstein-de Sitter universe has 
been given by a number of authors, including Matarrese et al (1998), (equation 
(6.10)), Bartolo et al (2005) (equation (8)), Hwang et al (2012), (equation (43)). 
We have shown that the expressions for (2)5 given by these authors can be written in 
our form with the arbitrary function C{x^) given by ( |39l) with g^n = | and a^i = 
in the first and third papers. 



^■'^See their equation (29). Note that their g{r]) equals our g{x) up to a constant multiple: 
g{x) = ^girj), and our f{x) equals their /{ij). Some rearrangement of the spatial dependence 
terms has to be done to relate their functions Ba{x^) to ours. 

^^See Tomita's equation (4.16). 
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3 Asymptotic behaviour at late times 

We now turn to the problem of determining the asymptotic behaviour of the second 
order perturbations as a; — )■ oo, that arise from a general (scalar) linear perturba- 
tion f l26ap . i.e. one that has both a growing mode ( and a decaying mode C7_. It 
is thus necessary to use the general solution fl26al) for We consider the three 
cases listed in the introduction. 



3.1 The de Sitter asymptotic regime 

In this section we determine the form of the nonlinear perturbation in the de Sitter 
asymptotic regime: x — oo with A > 0, i^' arbitrary (assuming background solutions 
that are forever expanding when i^' > 0). The first step is to use the asymptotic 
form of to determine the asymptotic form of the source term S{'^) in f l26cp . 
Then the asymptotic behaviour of can be determined using equation (I26ap . 



3.1.1 Linear perturbations 

The asjTuptotic form of is determined by the functions 'H(x) and I{x) through 
equation (^^. It follows from f lT3|) and (I19bp that as x — t- oo 



n = Xx{l + \kxx-^ + 0{x-'^)) , where A := JVLk^UI, (45a) 



X-'I{x) = /oo - \x-' + 0(a;"^), where I^ := A'^ j jjj^dx. (45b) 



oo 



1 



Equation fl24|) then leads to 

^(x, X*) = \mx ^^^G j^{x^)x'^ + 0(x"^), as x oo, (46a) 

where 

(iteA(xO := IU{x') + C_(xO. (46b) 
Here we have introduced the scaling parameters 

mx-= — = ^""'^ kx ■= -— = ^^'^ (47) 
A^ f2A,o X"^ fl\fi 

It follows from (1211) . fll5ap (146 ap and (1160 that the asymptotic behaviour as x — > oo 
of the matter gauge invariants is given by 

1)2 

«(5 = mA^('teA-3C + 0(x-^), (48a) 



7/(iV = -( C - k^x^GA) + 0{x-^). (48b) 
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3.1.2 Second order perturbations 

A particular solution for ^"^^ is given by fl26ap . In order to determine the asymptotic 
behaviour of ^'^^ as x — )■ oo we need to determine the asymptotic behaviour of the 
source term 5(x,a;*). The key result is that S{x,x'') = /(x*) + 0(a;~^). We now 
derive this result, obtaining an explicit expression for the leading order spatial term 
fix'). 

It follows from (I70al) and the asymptotic forms of '■^^ and T/'-^V (details in the 
appendix) that 

S = n^[^^ + AV{'^)]+0{x-^). (49) 

Since S{x,x^) is quadratic in and its derivatives it will have a multiplicative 
factor of m\. We thus rescale S{x,x^) and S(x, x*) by defining 

S{x, x') = 'IX^ml S{x, x'), S{x, x') = fA^m^ S{x, x'). (50) 

It now follows from fl46aD . fH9|) and ([50]) that 

S{x,x') = gA{x') + 0{x-^), (51) 

where 

6;A(xO:=|(('tei + 4I)(«GA)). (52) 

This in turn implies that 

E{x,x') = gAix')x + 0{l). (53) 

We now use fl26al) and ( |53l) to obtain the following asymptotic expansion for 
as X — J- oo (details in the appendix): 

(%(x, x') = fm^ ((2)G'a(xOx-i - gA{x')x-^ + 0{x~'')) , (54a) 

where 

This improper integral converges since ?/(x) = (9(x) and S(x, x*) = 0(x) as x ^ oo. 

The asymptotic behaviour as x — )■ cxd of the remaining second order perturbation 
variables is determined by equations ( IT5l) . (!70l) and (|T6l) : 

(2)$ = (2)^ + 0(x~2), (55a) 
(2)^ = m2^(2)^A - Smx^^^C + O(x-i), (55b) 
7^(2)v = -rriA (('^C - Ga) + 5^ + O(x-i), (55c) 



where 



5v = 25* 



r)2 



3C) D,(C-A;a«Ga) 



(55d) 



Here ^^X, as defined by ([90]), is a quantity that depends quadratically on the first 
order function ^^^GAix"^) given by (146bp . If the decaying mode is zero (C_(x*) = 0) 
then (2)^ depends quadratically on as does Sv In this special case one can also 
express '•^^Ga and ^'^X as a linear combination of the terms quadratic in ( and its 
spatial derivatives that appear in f l33]) and f l35]) . 
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3.2 The Milne asymptotic regime 

In this section we determine the form of the nonhnear perturbation in the asymptotic 
to Milne regime: x — > oo with A = 0, < 0. 

3.2.1 Linear perturbations 



It follows from (IT^ and (]19b|) that the asymptotic behaviour of H and I{x) as 
X — 7- oo is given by 

H = y/-K (l + ^rrikX^^ + 0(x~^)) , where := = — — , (56a) 

3 1 

(-K)2 Ux) = x{l- |mfc— + Oix-^)). (56b) 
Equation ( l24l) then leads to 

111 X 

\E'(a;,a;') = |mfcC(x')a;"-^(l - fmfc ) + 0(x"^), as x ^ oo. (57) 

It follows from fl2T]) . fll5ap . fl57|) and f|T6|) that the asymptotic behaviour as x — )■ oo 
of the matter gauge invariants is given by 

r)2 

= (mfe — -3)C + 0(x-Mnx), (58a) 
•H{i)v = -fm^C — + C'(x-^). (58b) 

3.2.2 Second order perturbations 

A particular solution for '•^Kl' is given by fl26ap . In order to determine the asymptotic 
behaviour of ^^•Kl' as x — )■ oo we need to determine the asymptotic behaviour of the 
source term iS(x, x*). The key result is that iS(x,x*) = /(x*)x~^ + 0(x~^lnx). We 
now derive this result, obtaining an explicit expression for the leading order spatial 
term /(x*). 

It follows from (l7Ual) and the asymptotic forms of \E' and "HV (details in the 
appendix) that 

S = ir[3^2 ^ 4p(^)] + 1[-(D^)2 + 4D2p(^)] + 0(x-3(lnx)2). (59) 

Since S'(x, x*) is quadratic in and its derivatives it will have a multiplicative 
factor of ml. We thus rescale S'(x,x*) and §(x, x*) by defining 

S{x, x') = 1{-K)ml S{x, x'), §(x, x') = 1{-K)ml S(x, x'). (60) 

It now follows from (I57D, (El) and (EOl) that 

S{x,x') = gk{x')x-^ + 0{x-%\nxf), (61) 

where 

Gkix') :=-l[3C + AViC) + ((DC)^ - AT^'V^Q)] . (62) 
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This in turn implies that 

S{x,x^)dx = hm x*), exists and is non-zero. 



(63) 



We now use the general solution (I26ap in conjunction with ( !63|) to obtain the 
following asymptotic expansion for ^'^^ as x — )■ oo (details in the appendix): 



(%(a;,x') 



3 2 —1 



In a; 



X 



+ 0{x~ 



where 



S{x, x^)dx. 



(64a) 



(64b) 



In the limit x — ?■ oo the remaining second order perturbation variables are de- 
termined in terms of (^)\E' and (^)\1> by equations ([15]), ([70]), ([28]) and f[T6l) : 



(2)$ = (% + 0(x 



-2n 



mkirrik 



-3)(2tefc + 0(x-Mn 



Inx 



X 



+ 0(x~^). 



(65a) 
(65b) 

(65c) 



We note that if the decaying mode is zero one can express as a linear com- 
bination of the terms quadratic in ( and its spatial derivatives that appear in f l33]) 
and ([35]). 



4 Discussion 

We have used the minimal system of governing equations for second order pertur- 
bations that we developed in UW3 to investigate the behaviour of nonlinear scalar 
perturbations of FL cosmologies with dust and cosmological constant as matter- 
energy content. Although we have given a new formulation of the solutions of 
the governing equations, our main contribution in this paper has been to analyze 
the properties of the solutions, and this has been facilitated by the new form of 
the solutions. In particular we have shown that there are significant differences in 
the asymptotic behaviour at late times in three cases, characterized by K and A, 
namely, the Einstein-de Sitter regime, the Milne regime and the de Sitter regime, as 
described in the introduction. These differences are best illustrated by considering 
the behaviour of the linear and nonlinear density perturbations ^^^6 and ^^•'<5. For 
ease of comparison we repeat the three asymptotic expressions for ^^^S as x — )■ oo: 

^'^SdeSitter = ^^^(/ooC + C^) - 3C + ©(x"!), (66a) 

^''^6 Milne = mfc — C - 3C + O (x-Mn x) , (66b) 
WdEdS = lx^( + Oil), (66c) 
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where I^o is given by f l45bp . Note that the expressions for the density perturbation 
^"^^S in the three asymptotic regimes are given by equations (]55bll . (I65bl) and ( l44l) . 

An arbitrary scalar perturbation depends on two spatial functions C{x'^) and 
C_(a;*), that correspond to the growing and decaying modes, and on the state of the 
background model at some initial time, which determines the following constants: 

2 i/2o ^m,0 7 ^fc,0 ^m,0 /cj-tN 

m = Hoiim.o, rnx = — — , kx = - — , m^. = — — . (67) 

"A,0 "A,0 "fc,0 



It follows from f l45bp that 



loo = Iooifnx,kx). (68) 

By inspection of equations fIBB]) we see the following dependencies of the asymptotic 
states for the first order density perturbation ^^^S: 



Einstein-de Sitter 
Milne 
de Sitter 



m, C, (69a) 
m,mkX, (69b) 
m,mx,kx,C,C- (69c) 



We see that ^^^d has the most comphcated structure in the de Sitter asymptotic 
regime, depending on the ratio of flm,o and flk,o to fl\fi and on the decaying mode 
C__. The situation is different as regards the second order density perturbation. 
Since ^'^^ and hence ^'^^S has an integral dependence on the source term S (see equa- 
tion f l29p ). the decaying mode enters into the leading term in all three asymptotic 
regimes. 
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A Analysis for second order perturbations 

A.l Source terms for second order perturbations 

In this section we give the full expressions for the source terms and then illustrate 
their structure in the special case when the decaying mode is zero. 

The source terms for a perfect fluid model with a purely scalar linear perturbation 
are given in UW3 (see equations (61)). After specializing to the case of dust and 
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rearranging the terms we obtaid^M 

S = {Hxd^^f + 4ir^2 - i ((D^)2 - A{T)Nf) 

+ 2 [n\xd^ + 3fiA + f^fc) + |D2] (2P(^) + ^I)(V)) , 

Sy = 2A-^n^ [^^ + 4r'(^)] + 2 [5^(7/VDi^) + 2D(7/V)] , 



(70a) 

(70b) 
(70c) 
(70d) 



Here and elsewhere we use the following notation. First, the second order spatial 
differential operators are defined by 



(71) 



where Dj denotes covariant differentiation with respect to the spatial metric 7,^. 
Second, we use the shorthand notation 

(DA)' := (D'=A)(DfcA), V{A) := 5^^(D,A)(D,A), (72) 

where A is a scalar field. Finally, we use the scalar mode extraction operators 

= B-^B\ S'^ = |D-2(d2 + SK)'^ (73) 

as defined in UW3 (see equation (85a)). 

A. 1.1 Zero decaying mode at linear order 

It follows from f l70ap that S{x, x*) can be written in the form 

Six, x') = T,{x)C + T2{x)V{Q + (T3(x)(DC)' + T,{x)B^V{Q) , 



where 



T, = {Ug)\U-^f-^^k). 

T2 = -mgf ((/ - 1)' - + ^k) + i^/' 

1 2 
3m 



= Wo' (1 + ^f- 

and / and g are defined in fl32|) . The remaining source terms are given by 



§9^ 



xm 



(4D\' - 5(DC)') + £ (/(/ - 1) - 4fi,) 



1 + ^^(0 



(74) 

(75a) 
(75b) 

(75c) 

(75d) 

(76a) 
(76b) 
(76c) 



^^In UW3 the source terms are not given a label 5,, but can be identified by comparing equa- 
tions (|14al) and with equations (61) in UW3, noting the difference in choice of variables as 
described in equation (|105p . We have also used the first equation in (jl5bp in obtaining these 
expressions. 
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A. 2 Asymptotic expressions for the source terms 
A. 2.1 Asymptotic to de Sitter 

In order to derive ( H9l) we have to examine f l70ap in detail. We begin by using 
( |T3|) . f l46ap and ( ]48bl) to conclude that Vi^^) and AViY) have the asymptotic form 



V{^) = A{x')x-^ + 0{x-^), AV{Y) = B{x')x-'^ + 0{x-^), (77) 

and 

n^xd^ + Sn^ + n^) + fD^ = AV(x9^ + 3) + 0(1). (78) 
It then follows that 

x'^ixd^ + 3)r'(^) = s2l?(^) + 0{x-^), x^{xd^ + 3)^I?(V) = 0{x-'^). (79) 



The desired result ^ follows when ([78]) and ([79]) are substituted in fl70al) J^1 The 
asymptotic forms as x — > oo of the other source terms given by (1701) are: 



= Qml [Gl - P(G'a)] + 0(x"3), (80a) 

5b = 0(x-^), (80b) 

Sy = tuxGax + AV{C- /caGa) + 0{x-^). (80c) 

The leading order term in So is quite complicated, and is not needed for our pur- 
poses. 

A. 2. 2 Asymptotic to Milne 

In order to derive ( 159]) we have to examine ( I70ap in detail. We begin by using 

([13]), ([57]) and f l58b]) to conclude that that P(^) and AV{Y) have the asymptotic 



form 

p($) = f{x')x-^ + 0{x-^), AViY) = g{x')x-\lnxf + 0{x-*), (81) 

and 

n^ixd, + + = {-K){xd, + 1) + + 0(x-2). (82) 
It then follows that 

{xd^ + l)r'(^) = -r'(^) + 0{x-^), {xda: + 1)AV{Y) = 0{x~\\nxy). (83) 



The desired result ([59]) follows when and are substituted in fl70aD J^1 The 
asymptotic forms as x — )■ oo of the other source terms given by ( 170]l . are: 



= 0(x-2), (84a) 
5^ = fmfc^m-^ [4(D^ + 3K)C^ - 5(DC)^] x'^ + 0{x"Hnx), (84b) 
Sy = Irrik [C + ^V{C)]x-^ + 0{x-\\nxf). (84c) 



i®The terms in the first line of (jTOa]) are 'H^^'^ + 0{x^^). 

^^The terms in the first hne of (|70a| are 3^:*^ - ^(D^')^ + 0{x^^ lux). 
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A. 2. 3 Einstein-de Sitter 

The exact expressions for tlie source terms given by (ITOll are 

5 = i[(D^)2 + 20D2r'(^)], (85a) 
5$ =4[^^- (85b) 
So = 2^-1[4D2^2 _ 5(D^)2], (85c) 

Sy = fV{^). (85d) 

A. 3 Asymptotic expansions for second order perturbations 
A. 3.1 Asymptotic to de Sitter 

We outline tlie derivation of equation fl54ap for It follows from f l26ap . f H7|) 

and (EOD that 

Since 

B{x,x') = xgAix') + 0{l), ^ = l + 0(x-2), (87) 

Ax 

(see equations ([13]) and ([53])) it follows that lim^^oo x^'^M/{x, x*) exists. With ^'^^Ga{x'-) 
defined as in f l54bp we obtain 



The integral on the right can be expanded using (IHTj) . leading to the desired equa- 
tion (I54ap . 

To calculate we need to extend the asymptotic expansion of S given by ( l53l) . 
We write 

^{x,x') = gK{x')x + ^Q + 0{x-^). (89) 

When (l89l) and fl80cp are substituted into (I28)) . after making use of ([50]) . we obtain 
the desired equation (I55cp . with 

(2)C:=So + 4P(C-fcAGA). (90) 
A. 3. 2 Asymptotic to Milne 

We outline the derivation of equation fl64ap for We begin by using ( l50l) to 

write fl26ap in the form 



x(%(x,xO = |m^^ / ^rfx, (91) 



M r s(x,x' 

where 

H{x) := {-K)-^'^H{x) = (1 + mfcX-^)i (92) 
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We next derive an asymptotic expansion for S(x, x*). It follows from equations (126bp . 
(EO]) and f lMbl) that 

POO 

Six,x') -^^^Gkix') = - S{x,x')dx. (93) 



Substituting the expansion fl6T]) of S{x,x'^) and evaluating the integral gives 

S{x,x') = '^^^Gkix') - gk{x')x-^ + 0{x-^lnx), (94) 
which in conjunction with fl92|) leads to 

= (^te, - (|m, + + 0{x-Hnx). (95) 

On substituting this expression in (|9T1) and evaluating the integral we obtain f l64al) . 

B The matter gauge invariants 

In this appendix we define the two types of matter gauge invariants that we use in 
this paper. We consider a perfect fluid with stress-energy tensor 

T% = {p + p)u''ub + pS%, (96) 

with a linear equation of state p = wp with w constant. We begin by reformulating 
the Replacement Principle for the stress-energy tensor of a perfect fluid, as given 
in UW3 (see equations (84)), by viewing it as a function of the variables F = 
{p,Va,gab), where Va = a~^Ua is the conformal fluid velocity and (jab = a'^Qah is 
the conformal metric tensor. The perturbations of the stress-energy tensor can be 
written symbolically in the form: 

(i>r«, = n((i)F), ^^"T% = n((2)F) + TM'-'^F), (97) 

where T% is the linear leading order operatoi@ and T% is the quadratic source term 
operator, and 

(^)F={M'^'^p^%a,^%b), (98) 
with r = 1, 2, denotes the perturbations of p, Va and (jab, while M. is defined by 

:=((0)p + (0)p)-V2. (99) 

For a linear equation of state the perturbations in the pressure are given by 

= w ("V, for r = 1, 2; (100) 

We associate gauge invariants with Ai^ ^'^^% and the variables '•^^F as follows. 
The gauge invariants *^''^y4[X] associated with the perturbations ^'^^A of an arbitrary 
tensor A are defined b\o 



(^U[X] := (^U - (101a) 
(2U[X] := - £c2)x^'^A - £a)x i^A^'^ " £wx^°'^A) . (101b) 



^°Note that the same operator T^f, acts on both '^^^F and ^^''F. 

^^See Nakamura (2007), equations (2.26)-(2.27), and UW3, equations (81). Here we are omitting 
the factor of a" in the latter equations. 
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In terms of this definition the gauge invariants associated with Ai^ ^^\r% and Ai^ ^"^^p 
are defined by replacing the tensor A in (11 Oil) by T% and by p: 

^'>F%[X] ■.= M^ '^''TWXl ^'^5[X] ■.= M'^ (102a) 

while the gauge invariants associated with '^'^ha = a~^^'^\ia and ^"^^fab '■= ^^^ab = 
c^"^ ^^''dab are defined by replacing the tensor A in fllOll) by Ua and by gab'- 

^%,[X] := a-' MwjX], ^Xb[X] := a'^ ^%ab[X]. (102b) 

The Replacement Principle states that the gauge invariants associated with A^^ ^'^'T% 
and with (''^F are related by the same operators as in (^71) : 

= n((^)F), = n(('>F) + n(^'^), (103a) 

where '•''^ is shorthand for 

(^>F[X] = {}%X]^%a[Xi '^XbiX]). (103b) 

Finally the scalar velocity perturbations are defined by 

^%[X] = S' ^%i[X], for r = 1, 2, (104) 

where the scalar mode extraction operator 5* was given in fl73l) . We can now state 
that the first set of matter gauge invariants for linear and second order scalar per- 
turbations are ^^^S[X] and '^'"V[X], with r = 1,2, relative to an arbitrary choice of 
gauge. 

Before continuing we briefly digress to comment on the relation between our 
variables and *^^)5[X] and '^^V[X] and the variables used by other authors. First, our 
V equals up to sign the variable v used by others@ Second, as follows from (p9|l 
and (I102ap . our 6 is related to the usual relative density perturbation S =*^''^p/*^°^p 
according to 5 = (1 + w)d, where w is defined by ^'^^p = w^'^^p. 

We also define matter gauge invariants directly in terms of the stress-energy 
tensor @ 

(^)V[X] := 5*(")T°[X], (106a) 
("to[X] := -S' (Di(")T°o[X] + 3H(")T° [X]) . (106b) 

In order to relate these gauge invariants to ^''^<5[X] and '^^V[X] we need to obtain 
explicit expressions for the components of '^^)T'6[X]. Referring to (!96|) and (!97|) we 

^^To make comparisons we restrict to the Poisson gauge {X — Xp) at linear order. Noh and 
Hwang (2004) have = -v[Xp] (see equation (297)), Nakamura (2007) has v'-^^ = v[Xp] (see 
equation (5.20)) and Malik and Wands (2009) have = v[Xp] (see equation (8.21)). 

^■^The matter gauge invariants defined in UW3, using the symbols A and V instead of D and V, 
are related to the ones defined here as follows: 

('■)A[X] ('■V[X] =yl('''V[X]. (105) 

For A see equations (45b) and (41c), and for V, see (45c). 
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perform an e-expansion for T%{e) , gab{^) , Ua{e) and p(e) and obtain 

Jl{(^^F) = -M' ('■V, n((')F) = -7^^%, (107a) 

7° = = (2(1 + w)M' (i)p + /oo) (107b) 

J\{^''-¥) = 3wM^ (''V, m^'^F) = Y'v,„ (107c) 

t.,(('-)F) = 0, %i^'^F) = v,„ (107d) 

where 

Vij := 2vi{vj - foj) , (107e) 
and r = 1,2 in the first column. We now apply the Replacement Principle obtain- 



T°o((^>F) = -M(5, n(('^) = -7^^v,„ (108a) 

7° ((-)F) = -('■V,, n((i>F) = (2(1 + w)d + foo) v„ (108b) 

J\{^'-^) = 3w r\{^'^) = (108c) 

t,,(('^>F) = 0, t,(('^)=v,,, (108d) 



where 



:=2vi(v, -fo,). (108e) 
When we specialize to linear perturbations that are purely scalar we havj^ 

^%i[X] = Di(%[X], Vi,[X] := 2Div[X] D,(v[X] - B[X]) (109) 

In this case it follows from fllOSp . fll06p and fll08p that the two types of gauge 
invariants are related as follows: 

^%[X] = (%[X], (110a) 

^%[X] = (%[X] - 2S' [((1 + w)S[X] - $[X])DiV[X]] , (110b) 

^%X] = + 3H(%[X], (110c) 

(2)(5[X] = + 3n^^\[X] - 2D*V[X]Di(V[X] - B[X]). (llOd) 

Finally in the Poisson gauge, which is given by B[Xp] = 0, equations ( IllOp . when 
specialized to dust simplify to give equations ( !T6ll . 
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